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We calculate AI = 3/2 kaon decay matrix elements using domain wall fermions and the DBW2 
gauge action at one coarse lattice spacing corresponding to a -1 = 1.3 GeV. We employ the Lel- 
louch and Liischer formula and its extention for non-zero total momentum to extract the infi- 
nite volume, center-of-mass frame decay amplitudes. The decay amplitudes obtained from the 
methods correspond to those from the indirect method with full order chiral perturbation theory. 
We confirm that the result is consistent with the previous result calculated with H-parity (anti- 
periodic) boundary condition by investigating the relative momentum dependence. We evaluate 
the decay amplitude ReA2 at the physical point by a chiral extrapolation with a polynomial func- 
tion of nv\ and the relative momentum as well as the AI = 3/2 electroweak penguin contributions 
to e'/e. We found that the result of ReA2 reasonably agrees with the experiment. 
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1. Introduction 

It is difficult to directly calculate the K — > %% weak matrix elements on lattice due to dif- 
ficulties of calculation of the two-pion state in finite volume. Indirect method [jl|] is one of the 
candidates to avoid the difficulty. In the indirect method K — > %% process is reduced to K — > % 
and K -> processes through chiral perturbation theory (ChPT). The RBC and CP-PACS [|] 
Collaborations calculated full non-leptonic kaon decay processes with the method. Their final re- 
sults of e'/e, however, have the opposite sign to the experiment. In the calculations there are many 
systematic errors, e.g., calculating at one finite lattice spacing with quenched approximation, and 
using the reduction with tree level ChPT. The indirect method might cause larger systematic errors 
than other sources, because the final state interaction of the two-pion is expected to play an im- 
portant role in the decay process. Thus we have to treat the scattering effect of the two-pion state 
directly on lattice to eliminate this systematic error. 

There are two main difficulties for the direct method, where the two-pion state is calculated 
on lattice. The one is to extract the two-pion state contribution with non-zero relative momentum 
in the K — > %% four-point function, which was pointed out by Maiani and Testa [HJ. So far there 
are several ideas for solving the problem. A method is employed with a proper projection of the 
K — ► 7T7T four-point functions [Q]. In the method we need complicated calculations and analyses, 
e.g., diagonalization of a matrix of the two-pion correlation functions [Q], to treat the two-pion 
state with non-zero relative momentum on lattice. A simpler idea, where complicated analyses 
are not required, is to prohibit the zero momentum two-pion state. Recently Kim [0] reported 
an exploratory study of the idea with H-parity (anti-periodic) boundary conditions in the spatial 
direction. He succeeded to extract the two-pion state with non-zero relative momentum from the 
ground state contribution of correlation functions, because the zero momentum two-pion state is 
prohibited by the boundary condition. We can also forbid the zero momentum two-pion state 
in center-of-mass (CM) frame by performing the calculation in non-zero total momentum (Lab) 
frame, \P\ ^ 0. In the frame the ground state of the two-pion is \tc(0)k(P)), which is related to the 
two-pion state with the non-zero relative momentum in the CM frame. Thus, we can extract the 
two-pion state with non-zero momentum from the ground state contributions [JSj] as well as in the 
H-parity boundary case. 

The other difficulty is related to the finite volume correction due to the two-pion interaction. 
We have to pay attention to the finite volume effect of the two-pion to obtain matrix elements in the 
infinite volume, because it is much larger than that of a one-particle state. Lellouch and Liischer 
(LL) [Q] suggested a solution which is a relation between the finite and infinite volume, center-of- 
mass frame decay amplitudes. However, this relation is valid only in the CM frame with periodic 
boundary condition in the spatial direction, so that we need a modified formula when we utilize 
H-parity boundary condition [Q] or Lab frame calculation. Recently, two groups, Kim et al. [10] 



and Christ et al. [|1 1|], suggested a formula which is an extension of the LL formula for the Lab 
frame calculation. 

Here we attempt to apply these two methods, Lab frame calculation and the extended for- 
mula, to the calculation of the Al = 3/2 kaon weak matrix elements. Our preliminary result was 
presented in Ref. [|T2]]. The results presented here are systematically larger than the preliminary 
results reported in last year lattice conference. During the past year we discovered that larger time 
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separations were needed to remove excited state contamination. 



2. Methods 



Through the extended formula [|R], 11], the infinite volume decay amplitude \A\ in the CM 
frame is given by 



\A\ 2 = Snf^J [p 



dS(p) d^p(q) 



M 



(2.1) 



dp dq 

where \M\ is the finite volume, Lab frame decay amplitude, 7 is a boost factor, E nn is the CM two- 
pion energy, and 8 is the scattering phase shift of the final state interaction. The relative momentum 
p is defined by the two-pion energy, p- = E 2 x /4 — m 2 n . The function (pp with P being the total 
momentum, derived by Rummukainen and Gottlieb []l3|], is defined by 



tan (j>p(q) 



3/2 



z p (i;<7 2 ;y) 



(2.2) 



where q 2 = (pL/2n) 2 , and 



2ob(i;r;y) 



4k J^ 3 n\ + n\ + r 2 («3 + 1/2) 2 - q 2 



(2.3) 



in the P = (0,0,27l/L) case. The formula eq. ( |2~I| ) is valid only for on-shell decay amplitude, i.e., 
E nn = rriK as in LL formula [g]. 

We calculate the four-point function for A/ = 3/2 K —> Tin decay with total momentum P = 
and (0,0,2n/L). The four-point function Gi{P;t,t n ,tK) is denned by 



Gi(P;t,t„,t K ) = mK^P-tK^O^^K+n-iP-t^O), 

3 /2 

where the operators O] are lattice operators entering AI = 3/2 weak decays 
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(2.5) 
(2.6) 



with (^)l = ^7^(1 — 75)^, (qq)r = ^7^(1 + 75)^» and a, ft being color indices. and are the 

operators in the (27,1) and (8,8) representations of SU (3)l® SU (3) r with 1 = 3/2, respectively. 

3/2 3/2 

O mg8 is identical to 8 g , except the color indices are mixed. We also calculate the four-point 
function of two pions G 7t7l (P;t,t7c) and the two-point function of the kaon Gx{P',t,tK) with zero 
and non-zero total momenta, to obtain the needed energies and amplitudes. 



l3/2. 



3. Simulation parameters 

We employ the domain wall fermion action with the domain wall height M = 1.8, the fifth di- 
mension length L s = 12 and the DBW2 gauge action with j8 = 0.87 corresponding to aT l = 1.31(4) 
GeV and m res = 1.25(3) x 10~ 3 . The lattice size is L 3 • T = 16 3 • 32, where the physical spa- 
tial size corresponds to about 2.4 fm. Our simulation is carried out at four u,d quark masses, 
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Figure 1: Measured values of scattering amplitude T. and 8 obtained from CM and Lab calculations. Solid 
lines are fit results. 



m u = 0.015,0.03,0.04 and 0.05, for the chiral extrapolation of the decay amplitudes, using 252 
gauge configurations, except for the lightest mass where we use 370 configurations. For interpo- 
lations of the amplitudes to the on-shell point, we calculate the decay amplitudes with six strange 
quark masses, m s = 0.12,0.18,0.24,0.28,0.35 and 0.44, at the heavier three m u , while we use the 
three lighter m s at the lightest m u . It is enough for the on-shell interpolation with the three m s in 
the lightest m u case. We will see the point in a later section. 

We fix the two-pion operator at t % = 0, and the kaon operator = 20 to avoid contaminations 
from excited states as much as possible. A quark propagator is calculated by averaging quark 
propagators with periodic and anti-periodic boundary conditions for the time direction to obtain 
a propagator with 2T periodicity. We calculate the quark propagators with Coulomb gauge-fixed 
wall and momentum sources. 



4. Results 

The final state of the AI = 3/2 kaon decay is the S-wave isospin 1 = 2 two-pion state. The 
scattering phase shift of the two-pion state can be obtained by the finite volume formulae for the 



CM 014Q and the Lab frames [ |13| ] with the two-pion energy in each frame. We define "scattering 
amplitude" T(p) = tan8(p)E K7[ /2p, where E %% is the two-pion energy in the CM frame. The 
scattering amplitude is used for the chiral extrapolation of the phase shift with a global fitting for 
m\ and p 1 with a polynomial function a^m^ + a2om A 7l + amp 2 J r a\\m 2 7l p 1 . While p 4 and m A % should 
both be treated as second-order in ChPT, for our calculation we have only two different relative 
momenta. For that reason we omit the additional p A term. The figure [T] shows the result of T(p) 
and the fit results at each pion mass in the left figure, and the measured values of 8(p) in the right 
figure. The phase shift at the physical pion mass m n = 0.14 GeV, plotted by solid line with the error 
band in the right graph, is comparable with the prediction from ChPT with experiment. In order to 
utilize the extension eq.(|2Tl|) of the LL formula, we evaluate the derivative of the phase shift from 
the fit result, while the derivative of the function 0p is obtained by a numerical derivative. 

All the off-shell decay amplitudes are determined by the ratio of correlation functions with 
i = 27,88 and m88, \Z3Gi{P;t,t n ,t K )Z %% Z K /G nw {P;t,t % )G K {P;t,tK), where Z n% and Z/f are the 
overlap of the relevant operator with each state. We determine the off-shell amplitude in the flat 
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Figure 2: Interpolation of off-shell decay amplitude to on-shell in CM and Lab frames. Closed(open) 
symbol is off-shell(on-shell) amplitude. 
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Figure 3: Infinite volume decay amplitude of 27 operator obtained from different calculation methods. 



region of the ratio as a function of t, because the ratio will be a constant for those values of t when 
these correlation functions are dominated by each ground state. 

The figure || shows the off-shell decay amplitude of the 27 operator in finite volume and its 
interpolation to the on-shell Ek = E n7t in the cases of both the frames. In the lightest pion mass 
case, we can linearly interpolate the data with the three kaon energies. On the other hand, the 
off-shell decay amplitudes at the heavier pion masses have large curvature for the kaon energy, so 
that we fit the data with a quadratic function. In all the pion mass cases, the interpolated, on-shell 
decay amplitude and the fit line are plotted in the figure. 

The decay amplitudes on finite volume in each frame are converted to those of the CM frame 
in the infinite volume through the LL formula []|] and its extended formula eq.(2J) using the 
derivatives of 8 and (f)p obtained in the above. In Fig. |3| we plot the infinite volume decay amplitude 
of the 27 operator obtained from the different frames at the lightest pion mass as a function of p 2 . 
The previous result obtained from H-parity boundary calculation is also plotted to compare with 
these results. The amplitude obtained from the Lab frame is consistent with the line determined 
from those of the CM frame and H-parity boundary calculations. 

The infinite volume decay amplitudes of all the operators are plotted in Fig. ^ The result of 
the 27 operator seems to vanish at the chiral limit with zero relative momentum, while the other 
elements remain a constant in the limit. These trends of the pion mass dependence are reasonably 
consistent with the prediction of ChPT at leading order [J8J. In order to investigate the nv\ and 
relative momentum p 2 dependences, we carry out a global fitting for each decay amplitude for m 2 n 
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Figure 4: Infinite volume decay amplitudes obtained at different relative momenta. Diamond symbols are 
fit results with m\—p 2 — 0. Open symbols are omitted in global fit. 
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Figure 5: Measured ReA2 and result from previous works. Open symbols are omitted in global fit. In left 
panel dashed lines denote fit results. 



and p 2 assuming a simple polynomial form as 

A o+A n) m 2 z +A ol p 2 +A n mlp 2 . (4.1) 

The # 2 /d.o.f. of the 27 operator is larger than 5 with all the four pion mass data, so that we 
exclude the heaviest pion mass data in all the operators from the following analysis to make the 
# 2 /d.o.f. reasonable. The constant term of the 27 operator is consistent with zero within the error, 
Aqq = —0.0020(29), as expected. In other operators, we employ the same fitting function as the 
27 operator except A u =0. The fit results are plotted in the figures, and those at the both limits, 
nv\ = p 2 = 0, are represented by diamond symbols. It should be noted that while the linear fit to 
\Axi | is the correct leading CHPT behavior, it is being used in a region of quite large mass. Likewise 
the linear fit to \A%% \ and |A m 88 1 omits possible logarithm terms which are of the same order. 

The weak matrix elements are determined using non-perturbative matching factors previously 
calculated with the regularization independent (RI) scheme in Ref. [0|. The electroweak con- 
tributions of e'/e, {Qi)i and {Q%)%, are evaluated from |Agg] and |A m gg | . At the physical pion 
mass m % = 0.14 GeV and relative momentum p = 0.206 GeV, we obtain the matrix elements 
(<2?>2(a0 = 0.2473(64) and (gghG") = 1-160(31) GeV 3 at the scale n = 1.44 GeV using the 
same fitting form as in the |Agg| and \A m gg\ cases. 
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We calculate ReA2 from the weak matrix elements in the RI scheme with the Wilson coef- 
ficients evaluated by NDR scheme calculated in Ref. [0]. The left figure in Fig. || shows ReA2 
obtained from the CM and Lab frame calculations. The measured values of ReA2 are almost same 
as those of the 27 decay amplitude, in Fig. |], apart from the overall constant. The reason is that 
the main contribution of ReA2 comes from the 27 amplitude. We carry out a global fit using the 



same polynomial assumption eq.(|4.1|) except the constant term A 00 = to evaluate the result at the 
physical pion mass m n = 0.14 GeV and momentum p = 0.206 GeV. In the fitting we omit the data 
at the heaviest pion mass because of the same reason mentioned in the above. We can also carry 
out a reasonable fitting, # 2 /d.o.f. = 1.2, without An. The fit results at each pion mass and the 
physical point are plotted in the figure. 

We plot the results of ReA2 at the physical point in the right panel of Fig. || together with those 
of the previous works using the indirect method [§, [|] and direct calculation with evaluating the 
finite volume, two-pion interaction effect through ChPT [|l5|]. We estimate ReA2 = 2.26(41) and 
1.61(24)x 10~ 8 GeV at the physical point with and without An term, respectively. These results 
agree with each other and the experiment within two standard deviations. The result is encouraging, 
albeit it includes systematic errors due to quenched approximation, finite lattice spacing effects, and 
heavy pion and kaon masses. 
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